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of  multiple  reflections  between  two  simple  bodies,  two  adjacent  parallel  con¬ 
ducting  cylinders,  was  undertaken  to  determine  the  effect  of  multiple  reflec-  * 
tions  on  the  image  and  the  dependence  on  polarization  of  the  incident  field 
and  the  receiving  antenna  or  antennas  used  to  record  the  scattered  field. 

Cylinders  are  chosen  to  make  the  problem  analytically  tractable  and  to  provide 
analytical  results  against  which  experimentally  obtained  results  can  be  compared. 
Excellent  Agreement  between  experimental  and  analytical- results  was  obtained. 

The  results  show  that  cross-polarized  measurement,  when  the  object  is  illuminated 
for  example,  by  right-hand  circularly  polarized  (RHCP)  plane  wave  and  the  receiving 
antenna  measuring  the  scattered  field  is  left-hand  circularly  (LHCP) ,  can  be  used 
to  suppress  the  effect  of  multiple  reflections  between  the  cylinders  and  to  eli¬ 
minate  noticable  artefact  in  the  reconstructed  image.  This  result  can  be  readily 
interpreted  by  noting  that  reflections  are  specular  and  that  the  state  of  polari¬ 
zation  is  reversed  after  each  specular  reflection.  If  the  receiver 
polarization  on  the  other  hand  is  RHCP,  artefacts  due  to  multiple  scattering  be¬ 
tween  the  cylinders  will  appear  in  the  image.  The  artefact  does  not  correspond 
to  any  physical  detail  on  the  object,  however,  it  can  convey  additional  informa¬ 
tion  about  the  imaging  geometry.  This  study  also  illustrates  well  the  power  of 
microwave  diversity  imaging  as  a  tool  for  the  understanding  of  electromagnetic 
scattering. 

A  study  of  microwave  diversity  imaging  in  the  presence  of  severe  clutter  as 
applicable  to  space-borne  imaging  of  low  flying  objects  and  imaging  of  penetrable 
objects  has  also  been  undertaken.  The  results  show  that  excellent  image  dis¬ 
crimination  between  the  object  and  ground  clutter  can  be  realized  provided  that 
sufficiently  wide  spectral  and  angular  windows  are  utilized. 

We  have  also  found  that  microwave  diversity  imaging  furnishes  images  that  have 
attributes  that  make  them  ideal  for  use  in  automated  target  recognition.  These 
attributes  include:  (a)  the  image  is  always  centered  in  the  image  plane,  (b)  the 
image  is  edge  enhanced  and  hence  resembles  a  "primal  sketch"  where  detail,  important 
for  characterizing  and  recognizing  the  object  is  retained  and  all  other  detail, 
such  as  flat  surfaces,  is  discarded  reducing  thereby  the  dimensionality  of  the 
signals  being  dealt  with  which  is  highly  desirable  for  automated  recognition,  (c) 
although  image  information  (Fourier  space  data  )  in  microwave  diversity  imaging 
is  obtained  when  looking  at  the  object,  for  example  from  below,  at  a  slant  angle, 
the  displayed  image  detail  appears  as  if  the  object  is  observed  directly  from 
below.  This  is  a  consequence  of  the  tomographic  nature  of  microwave  diversity 
imaging.  All  of  the  above  features  make  the  image  information  obtained  ideal  for 
the  generation  of  target  classifiers  that  are  suitable  for  use  in  automated 
recognition  schemes  that  might  utilize  associative  memories  based  on  models  of 
neural  networks.  For  this  reason  we  have  initiated  a  study  of  optical  implemen¬ 
tations  of  neural  net  models,  where  robustness  and  fault  tolerant  features  of 
associative  recall  in  these  networks  are  combined  with  the  parallelism  and 
massive  interconnection  attributes  of  optics  to  provide  a  new  approach  to  optical 
processing  of  sensory  data. 
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RESEARCH  IN  IMAGE  UNDERSTANDING  AS  APPLIED  TO 
3-D  MICROWAVE  TOMOGRAPHIC  AND  PROJECTIVE  IMAGING  WITH  NEAR  OPTICAL  RESOLUTION 


1.  INTRODUCTION 

. —  4 

^  The  objective  of  research  under  this  grant  is  to  achieve  fundamental 
understanding  of  the  dyyamics  of  coherent  and  incoherent  image  formation  and 
inverse  scattering  employing  spectral,  angular  (aspect  dependent)  and  polari¬ 
zation  diversity  as  means  for  realizing  cost-effective  imaging  of  3-D  objects 
with  near  optical  resolution.  Special  attention  is  to  be  given  to  3-D  micro- 
wave  diversity  radar  imaging.  The  approach  is  Theoretical  -  through  develop¬ 
ment  of  a  unified  theory  of  microwave  diversity  Imaging  for  both  stationary 
and  non-stationary  objects  which  has  been  reported  on  extensively  in  preceeding 
reports  and.  Experimental  -  through  development  and  verification  of  methodol¬ 
ogies  for  efficient  microwave  diversity  imaging  that  can  be  employed  in  future 
broad-band  imaging  radar  networks  capable  of  providing  3-D  projective  or 

tomographic  images  of  remote  aerospace  targets  with  near  optical  resolution.  — 

\ 

To  date  unprecedented  image  resolution  has  been  achieved  (~1.5  cm)  on 
scale  models  of  relevant  targets  employing  (6-17)  GHz  spectral  window,  polariza¬ 
tion  diversity,  and  (~ir/2)  angular  aperture  using  a  specially  developed  experi¬ 
mental  microwave  imaging  facility  for  data  acquisition  and  either  Fourier 
inversion  or  an  equivalent  filtered  back-projection  reconstruction  algorithm  for 
image  reconstruction.  We  expect  the  results  of  this  effort  to  tell  us  how  to 
achieve  centimeter  resolution  on  remote  aerospace  objects  cost-effectively  using 
microwave  diversity  radar  networks  that  will  also  be  ultimately  capable  of  high 
speed  automated  target  recognition  and  classification. 


During  the  period  of  this  report  our  research  effort  focused  on:  (a)  the 
study  of  polarization  and  multiple  scattering  effects  in  microwave  diversity 
imaging,  (b)  imaging  in  the  presence  of  severe  clutter,  and  (c)  neural  models 
for  fault  toleranct  computing  and  robust  signal  processing  and  their  optical 
implementation.  Results  of  this  research  are  summarized  in  Section  2  of 
this  report  with  details  being  included  in  Section  6  as  appendices.  Discussion 
of  research  findings  and  of  new  research  directions  is  given  in  Section  3. 

Lists  of  publications,  theses,  and  other  activities  during  the  report  period 
are  given  in  Section  4. 

2 .  CURRENT  PROGRESS 

To  date  our  efforts  have  focused  on  the  investigation  of  methods  that 
enable  the  formation  of  highest  resolution  microwave  and  millimeter  wave 
Images  of  scattering  bodies  cost-effectively.  At  such  high  resolutions, 
the  effect  of  multiple  scattering  in  images  of  complex  shaped  scattering 
objects  become  particularly  important  as  an  image  degrading  mechanism  or  as 
an  effect  that  its  proper  understanding  and  interpretation  can  convey 
additional  information  about  the  object.  The  effect  of  multiple  scattering 
is  to  impose  on  the  microwave  image  artefact  or  detail  that  makes  it  depart 
from  the  visually  expected  or  "optical"  image.  For  this  reason,  a  systematic 
theoretical  and  experimental  study  of  multiple  scattering  effects  in  the 
imaging  of  simple  bodies  e.g.,  consisting  of  two  conducting  or  dielectric 
cylinders  and  the  role  played  by  polarization  in  such  effects  was  conducted 
(See  Appendix  I  for  details) .  Analysis  of  multiple  scattering  effects 
between  cylinders  explained  well  the  artefacts  observed  in  earlier  imaging 
experiments  of  such  objects  and  their  dependence  on  polarization  of  the 


illuminating  and  receiving  antennas.  The  advantage  of  using  circularly  polarized 
illumination  and  cross-polarized  illuminating  and  receiving  antennas  to  suppress 
such  artefacts  is  demonstrated.  This  work  also  illustrates  dramatically  the  power 
of  microwave  diversity  imaging  techniques  as  a  tool  in  study  and  understanding  of 
electromagnetic  scattering. 

Another  task  studied  during  the  period  of  this  report  is  evaluation  of 
the  capabilities  of  microwave  diversity  imaging  in  the  formation  of 
recognizable  images  in  the  presence  of  severe  background  clutter.  This 
problem  is  Important  for  certain  modalities  of  radar  imaging  such  as  space- 
borne  Imaging  of  low-flying  aircraft  employing  a  "spot -light"  imaging  mode 
and  imaging  of  internal  structure  of  penetrable  dielectric  bodies  in  the 
presence  of  strong  surface  reflections.  Details  and  results  of  this  study 
are  given  in  Appendix  II.  There  a  projection  image  of  the  visible  scatter¬ 
ing  centers  of  a  test  object,  a  72:1  metalized  scale  model  of  the  space 
shuttle,  was  obtained  employing  a  spectral  window  of  11  GHz  (6.1-17.5  GHz) 
and  an  angular  window  of  90®.  The  shuttle  model  was  imaged  when  situated 
40  cm  above  corrugated  aluminum  ground  planes  of  different  roughness  to 
simulate  ground  clutter  from  different  terrain.  Excellent  polarization 
enhanced  Images  of  the  target  shape,  clearly  delineated  above  the  image 
of  the  ground  plane  were  obtained.  These  compare  very  well  with  flood-lit 
side-view  photographs  of  the  configuration  showing  the  silhouette  of  the 
object.  The  results  clearly  demonstrate  the  advantages  of  using  extended 
spectral  and  angular  windows  in  this  mode  of  imaging. 

Another  research  task  considered  during  this  period  was  the  study 
of  neural  models  and  their  optical  implementation  for  achieving  associative 
memories  that  can  be  employed  in  robust  fault -tolerant  computing  and 
Information  processing  of  sensory  (radar,  optical,  infra-red,  etc.)  data. 

The  results  of  optical  implementation  of  a  32  neuron  network  and  its 
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evaluation  as  an  associative  memory  are  detailed  in  Appendix  III.  An 
outer  product  storage  algorithm,  the  so  called  Hopfield  model,  is  employed 
to  realize  distributed  information  storage  akin  to  holographic  storage 
of  information.  This  work  demonstrates  clearly  the  power  of  optics, 
namely  parallelism  and  ease  of  making  global  or  semi-global  interconnection 
in  the  implementation  of  simplified  models  of  neural  nets.  On  the  other 
hand  neural  net  modeling  introduces  the  power  of  nonlinear  feedback  into 
optics  to  result  in  optical  processors  that  can  process  sensory  information 
collectively  with  a  degree  of  robustness  and  fault  tolerance  that  was  not 
possible  earlier. 

3.  DISCUSSION  AND  NEW  RESEARCH  DIRECTIONS 

The  principles  and  methodologies  of  microwave  diversity  imaging 
studied  in  our  work  show  that  it  is  possible  to  produce  high  resolution 
projection  or  tomographic  Images  of  radar  targets  with  near  optical 
resolution.  Such  Images  are  perfectly  suited  for  recognition  by  the 
eye-brain  system.  They  have  however  several  features  that  make  them 
ideal  for  use  in  automated  recognition  schemes.  These  include:  (a)  the 
reconstructed  image  is  always  centered  about  the  origin  of  the  image 
plane.  This  is  a  direct  consequence  of  the  target  derived  reference 
method  employed  to  produce  synthetically  during  pre-processing  of  the 
collected  data  a  phase  reference  point  centered  on  the  target,  (b)  the 
Image  is  edge  enhanced  containing  information  about  edges,  protrusions 
and  other  characteristic  detail  necessary  for  recognition  while  all 
other  detail  such  as  flat  surfaces  that  is  not  necessary  for  recognition 
is  discarded.  This  represents  a  highly  desirable  reduction  of  the 
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amplitude  and  plotting  the  sinogram  over  0  <  <f>  <  2tt  .  The  data  in  the 
sinogram  is  exactly  the  same  as  the  data  in  the  image  since  one  way  of 
forming  an  image  is  by  filtered  back-projection  [1]  where  the  range 
profiles  are  plotted  in  polar  format  and  backproj ected  to  form  the 
image  by  coherent  addition  of  the  contribution  of  the  various  back- 
projections  at  each  point  in  the  image  plane.  We  expect  sinogram 
classifiers  to  be  of  use  as  training  sets  for  two-dimensional  versions 
of  the  one-dimensional  associative  memory  described  in  Appendix  III. 

Such  memories  can  be  viewed  then  as  "filters"  that  respond  only  if  one 
of  the  entities  stored  in  them  or  a  close  version  of  one  of  them  (close 
in  the  Hamming  sense)  is  used  to  stimulate  or  initialize  the  memory. 

In  this  fashion  we  would  realize  a  rudimentary  automated  recognition 
system  that  is  both  robust  and  fault-tolerant  providing  thereby  a  new 
approach  to  automated  radar  target  recognition. 
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APPENDIX  I 


POLARIZATION  EFFECTS  IN  MICROWAVE  DIVERSITY  IMAGING 
OF  PERFECTLY  CONDUCTING  CYLINDERS 


Tah-Hsiung  Chu 

RCA  David  Samoff  Research 

Center 

Rm  3-231,  P.O.  Box  432 
Princeton,  NJ  08540 
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ABSTRACT 


A  theoretical  vector  development  of  the  directly  and  mutually 
scattered  wavefield  of  two  cylinders  in  a  bistatic  measuring  system  is 
presented.  The  two  scattered  fields  are  shown  to  be  polarization  depen¬ 
dent,  the  images  reconstructed  from  the  directly  or  mutually  scattered  far 
field  of  each  polarization  include  different  features  of  the  scattering  ob¬ 
ject.  An  automated  microwave  imaging  system  employing  frequency,  polari¬ 
zation,  and  angular  diversity  is  utilized  to  verify  the  theoretical  results . 
Analytical  and  experimental  results  sure  shown  to  be  in  good  agreement.  The 
polarization  effects  of  multiscatterers  are  also  discussed.  The  results 
presented  illustrate  the  utility  of  imaging  methods  in  the  study  of  electro 
magnetic  scattering  problems  and  in  image  understanding  in  general. 
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I.  INTRODUCTION 

The  use  of  polarization  diversity  over  an  extremely  broad  spectrum 
in  coherent  scattering  measurements  as  a  means  of  acquiring  more  informa¬ 
tion  about  the  scattering  object  has  been  the  subject  of  extensive  study 
at  the  Electro-Optics  and  Microwave-Optics  Laboratory  at  the  University  of 
Pennsylvania  [1-4].  An  excellent  review  of  the  role  of  polarization  in 
conveying  useful  information  about  the  scatterer's  characteristics  is 
found  in  [5],  and  a  discussion  of  the  polarization  effects  and  symmetry 
is  given  in  [6].  In  this  paper  the  results  of  a  detailed 
examination  of  polarization  effects  associated  with  the  far  field  scat¬ 
tered  from  two  conducting  cylinders  is  presented.  The  influence  of  polari¬ 
zation  on  microwave  images  of  the  cylinders  and  on  the  occurrence  of 
image  artefact  is  determined  analytically  and  verified  experimentally. 

A  theoretical  vector  development  of  the  far  field  scattered  by  per¬ 
fectly  conducting  cylinders  in  a  bistatic  measuring  system  is  derived  in 
Section  II.  The  approach  is  first  to  derive  an  expression  for  the  field 
scattered  from  a  single  cylinder  illuminated  by  a  right-hand  circularly 
polarized  (RHCP)  plane  wave.  This  expression  is  then  extended  in  Section 
III  to  the  case  of  the  directly  scattered  (or  first  order)  field  of  two 
cylinders.  The  mutually  scattered  (or  second  order)  field  from  two  cylin¬ 
ders  is  then  formulated  in  terms  of  a  coherent  summation  of  the  scattered 
field  from  each  cylinder  due  to  the  incident  wave  scattered  from  the  other 
cylinder.  Theoretically,  one  should  proceed  in  this  fashion  to  calculate 
the  n-th  order  scattered  field.  However  those  higher  order  scattering 
terms  are  found  to  become  negligibly  small  when  the  spacing  between  the  two 


cylinders  is  large  compared  to  their  radius  and  to  the  illuminating 
wavelength.  Therefore  in  order  to  simplify  computational  tasks .only  the 
second  order  scattered  field  is  analytically  derived  in  this  paper.  This 
is  called  the  mutually  scattered  field.  The  total  wavefield  scattered  by 
two  cylinders  becomes  the  coherent  superposition  of  the  directly  scattered 
component  and  the  mutually  scattered  component  for  each  polarization.  In 
Section  IV,  the  use  of  an  automated  microwave  imaging  facility  employing 
frequency,  polarization  and  angular  diversity  to  verify  the  results  of 
theoretical  analysis  is  described.  The  analytical  and  experimental  results 
are  shown  to  be  in  good  agreement.  The  results  show  that  depolariza¬ 
tion  effects  in  the  mutually  scattered  field  component  are  more  pronounced 
than  in  the  directly  scattered  component.  Depolarization  is  taken  to  be 
the  term  describing  the  change  of  the  polarization  of  an  electromagnetic 
wave  from  one  state  to  another  brought  about  by  the  scattering  process. 

The  notation  used  in  describing  the  state  of  polarization  in  this 

paper  is  given  next.  Referring  to  Fig.  1,  consider  a  plane  containing  the 

z-axis  perpendicular  to  the  plane  of  the  paper  and  the  incident  wave 
A  i  i, 

vector  k^.  The  quantities  E^  and  E^  are  orthogonal  components  of  a 

linearly  polarized  incident  plane  wave  lying  respectively  parallel  and 

perpendicular  to  this  plane.  When  right-hand  and  left-hand  circularly 

polarized  fields  E^  (RHCP)  and  E^(LHCP)  are  used  they  can  be  related  to 

E*  and  E*  by  [7] , 
z  <p 


Similarly,  we  can  express  the  right-hand  E*  and  left-hand  E®  circularly 


polarized  components  of  the  scattered  field  in  terms  of  the  scattered 
s  s 

fields  E  and  E ,  as 
z  <p 


F1  oi 

where  the  additional  matrix  accounts  for  the  coordinate  system 

L°  -iJ 

changing  "polarization  handedness"  upon  reflection  from  the  scatterer. 
Furthermore,  because  of  the  change  in  "handedness",  E?  and  Ej*  will  be 
called  the  co-polarized  and  cross -polarized  (or  depolarized)  components 
of  the  scattered  wavefield  when  the  incident  wave  is  RHCP  namely  E*. 

II.  POLARIZATION  EFFECTS  IN  SCATTERING  FROM  A  SINGLE  INFINITELY  LONG 
CYLINDER 

A.  Formulation 

Consider  an  infinitely  long  cylinder  of  radius  a  illuminated  by  a 
RHCP  plane  wave  E*  as  depicted  in  Fig.  1.  The  cylinder  axis  is  in  the  z 

K 

direction.  From  eq.  (1)  the  incident  electric  field  can  be  expressed  as 

E»  -  £  «}  -  K'  <3) 


where  from  [7,8] 
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are  the  parallel  and  perpendicular  polarized  components  of  the  incident 
wave,  Eq  *  n0Ho'T1o  t^le  intrinsic  impedance  of  free  space,  ^  and  (j)^ 

are  the  angles  between  the  x-axis  and  the  directions  of  incidence  and 
observation,  and  4>r  «  in  a  bistatic  measurement  system. 


Equations  (4)  and  (5)  can  be  expressed  in  a  Fourier-Bessel  series  as 


[7.  8], 
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where  J  (x)  is  the  Bessel  function  of  the  first  kind  of  order  n,  and  E 
n  o 

and  Hq  are  the  electric  and  magnetic  field  amplitudes  of  the  incident 
plane  wave. 

s  s 

The  scattered  wavefield  components  E^  and  E<j>  determined  from  the 
boundary  conditions. 
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where  R  designates  radial  distance  (see  Fig.  1),  are  given  by. 
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where  the  sunmations  over  n  here  and  in  subsequent  equations  is  understood 

(2) 

to  range  £rom  -«  to  <*>  ,  H  (x)  is  the  Hankel  function  of  the  second  kind, 

n 

(2) 

and  H*  (x)  is  its  derivative  with  respect  to  the  argument  x. 
n 

If  kR»l,  i.e.,  under  the  far  field  approximation,  the  bistatic  scat¬ 
tered  field  ((Jk  ?  0,  and  ^  «  IT  —  becomes 
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Using  the  transform  relationship  given  in  eq.  (2) ,  the  scattered  fields 
expressed  in  terms  of  circular  polarization  are 
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Equations  (15)  and  (16)  indicate  that  the  far  field  scattered  from  a 
cylinder  illuminated  by  a  RHCP  plane  wave  is  elliptically  polarized  in  the 
bistatic  case  and  also  in  the  monostatic  scattering  case  (<p^=0  and  4>r=T)  be- 
s  s 

cause  and  in  eqs.  (13)  and  (14)  differ  in  magnitude  and  are  not  in 
phase.  The  observed  scattered  far  field,  measured  for  example  by  a 
coherent  receiver,  can  be  separated  into  two  parts:  one  containing  the 
range  term  exp(-jkR)//R  appearing  outside  the  summation  terms  which  accounts 
for  the  propagation  time  delay  or  phase  shift  from  the  coordinate  origin 
to  the  receiver;  and  the  summation  terms  which  can  be  viewed  as  the  equi¬ 
valent  to  the  "Bojarski  identity"  [9,10]  which  strictly  speaking  conveys 
information  about  the  object  in  the  high  frequency  limit  (ka»l)  . 

It  was  shown  by  Bojarski  and  Lewis  assuming  physical  optics  and  Born 
approximations  hold  and  neglecting  polarization  effects,  that  a  Fourier 
transform  relationship  exists  between  the  shape  of  a  perfectly  conducting 
object  and  its  range-corrected  backscattered  far  field  frequency  response. 

In  the  next  section,  numerical  examples  are  given  to  examine  and  elucidate 
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the  Fourier  inversion  reconstruction  for  each  of  the  circularly  polarized 
fields  given  in  eqs.  (15)  and  (16)  in  the  microwave  regime.  Note  that  in 
arriving  at  the  results  in  eqs.  (15)  and  (16)  the  physical  optics  approxi¬ 
mation  was  not  involved  and  hence  they  hold  for  any  value  of  the  wavenumber  k 


Since  the  range  resolution  is  inversely  proportional  to  the  spectral 
range  of  the  illuminating  signal,  the  ability  to  pinpoint  the  reflection 
location  on  the  scattering  cylinder  is  greatly  enhanced  by  gathering  data 
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over  a  broad  frequency  range.  A  long  cylinder  can  therefore  be  used  as  a 

calibration  object  for  characterizing  broadband  microwave  scattering  system 

used  in  frequency  and  polarization  diversity  measurement  [4] .  In  Section 

IV,  the  scattered  field  E*  and  e5  from  two  cylinders  will  be  measured  in 

R  L 

an  anechoic  chamber  using  a  single  cylinder  as  the  calibration  object  for 
system  characterization. 

Ill .  POLARIZATION  EFFECTS  IN  SCATTERING  FROM  TWO  INFINITELY  LONG  CYLINDERS 
A.  Formulation  of  Direct  Scattering 

The  directly  scattered  wavefield  from  two  conducting  infinitely  long 
parallel  cylinders  illuminated  by  a  RHCP  plane  wave  as  shown  in  Fig.  3  can 
be  obtained  by  superposition  using  eq3.  (11)  and  (12).  This  leads  to  a 
summing  of  eqs.  (11)  and  (12)  with  each  multiplied  by  a  phase  term  due  to 
the  offset  distance  b  from  the  coordinate  origin  to  each  cylinder.  The 
resultant  field  expressions  referring  to  the  quantities  shown  in  Fig.  3 
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where  the  superscript  1  is  used  here  to  designate  the  direct  scattering, 
as  the  superscript  2  will  be  used  later  to  designate  the  doubly  scattered 
field  or  mutual  scattering  case.  The  subscripts  +  and  -  are  used  to  dis¬ 
tinguish  the  two  cylinders  and  0__^  is  the  angle  between  the  line  connect¬ 
ing  the  centers  of  the  two  cylinders  and  the  x-axis. 


In  the  far  field  approximation,  kR»l,  we  can  use  the  following  approxi 


(nations 
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Therefore,  after  normalizing  the  range  term  exp (- jkR) //r,  eqs.  (17)  and 
(18)  become 
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s  s 

where  ER  and  E^  are  the  expressions  for  the  scattered  fields  from  one 
cylinder  as  given  by  eqs.  (15)  and  (16)  excluding  the  range  terms.  Eqs. 
(23)  and  (24)  are  equal  to  eqs.  (15)  and  (16)  (excluding  the  range  term) 
multiplied  by  a  sum  of  two  linear  phase  terms  contributed  from  the  pro¬ 
jected  distances  +b  cosS^  coscfK  from  each  cylinder  to  the  coordinate 
origin . 


B.  Numerical  Results 

Figure  4  shows  the  range  profiles  derived  from  the  LHCP  and  the  RHCP 
components  of  the  directly  scattered  field  for  the  case  of  8  f  =  0°  and 
a*3.8cm,  b»*13.7cm.  These  are  near  replicas  of  single  cylinder  profiles 
displaced  by  distance  +b  cos4k  . 

By  measuring  the  backscatterad  field  components  for  different  closely 
spaced  values  of  8  . ,  correcting  them  for  the  range  factor  exp(-jkR)//R, 

-T 

and  presenting  them  in  polar  format  in  which  the  polar  angle  is  0  and 
radial  distance  represents  frequency,  we  can  obtain  a  slice  in  the  Fourier 
space  of  the  object  [3,4].  Fourier  inversion  of  the  slice  data  after 
interpolation  onto  a  rectangular  format  using  a  four  nearest  neighbor 
interpolation  algorithm  [12]  yields  in  accordance  to  the  projection-slice 
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theorem  [13)  a  projection  image  representing  the  projection  of  the  scat¬ 
tering  centers  of  the  object  on  a  plane  parallel  to  the  plane  of  the 
Fourier  space  slice. 
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Figure  5 (a)  and  <b)  shows  Fourier  plane  slices  obtained  from  the  E 

L 

1  s 

and  E  components  presented  in  a  polar  format  consisting  of  256  radial 
R 

lines  equally  spaced  over  360° ,  with  each  line  containing  64  equally 

spaced  frequency  points  covering  (6-17)  GHz  range.  Considering  that  the 

radial  (frequency)  and  angular  (aspect)  sampling  criterion  needed  to  avoid 

aliasing  are  respectively  6f  2.  c/4b  and  50  Amin/2b  [4],  the  acquired  data 

in  the  above  polar  format  is  somewhat  oversampled.  The  reconstructed  images 

for  each  circularly  polarized  component  as  obtained  via  a  two-dimensional 

Fourier  transformation,  are  given  in  Fig.  6.  Several  observations  can  be 

X  s 

made  based  on  the  images  shown.  First,  the  E  component  gives  two  dis- 

L 

tinct  circles  representing  the  external  cross-sectional  outline  images  of 

the  two  cylinders.  The  radii  and  spacings  of  the  cylinders  in  the  two  images 

correspond  quite  closely  to  the  values  a  and  b  of  the  object.  This  result 

agrees  with  the  physical  optics  diffraction  theory  given  by  Bo jar ski  and  Lewis 

X  s 

in  the  context  of  their  inverse  scattering  theory  [9,10].  However  the  ER 
component,  ignored  in  the  scalarized  physical  optics  approximation,  also  gives 
two  faint  circles  which  are  contributed  from  the  weak  depolarization  and 
creeping  wave  effects.  In  the  next  section,  we  will  analyze  the  mutual 
scattering  between  two  cylinders  in  each  polarization  state. 

C.  Formulation  of  Mutual  Scattering 

First  we  consider  the  case  when  the  cylinder  is  illuminated  by 
the  directly  scattered  wave  from  cylinder  to  obtain  the  contribution 


of  the  cylinder  to  the  second  order  mutually  scattered  wavefield.  The 


contribution  of  the  "+"  cylinder  to  the  mutually  scattered  wave  by  re-scat- 
tering  the  directly  scattered  wave  from  cylinder  can  be  obtained  by 
interchanging  the  "+"  and  subscripts  in  the  resulting  expressions. 

The  scattered  wave  from  the  cylinder  incident  on  the  cylinder 
is  rewritten  from  the  first  terms  in  eqs.  (17),  (18)  and  (5)  as 
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where  R+  in  accordance  to  Fig.  3  is  now  the  distance  between  the  center  of 
the  "+"  cylinder  and  an  observation  point  on  the  surface  of  the  second 
cylinder.  The  second  order  mutually  scattered  wavefield  from  the 
cylinder  can  then  be  expressed  in  terms  of  the  cylindrical  wave  function 
[7,8]  as 
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are  found  where  *A  and  *B  are  given  in  eqs.  (27)  and  (28).  Note  eq. 
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(34)  which  is  the  scattering  coefficient  of  the  _E^  component  is  the  same 

as  that  derived  by  Twersky  [15]  for  the  second  order  scalar  scattering 

2 

coefficient  of  two  cylinders.  The  other  scattering  coefficient  in 

2  s 

eq.  (35)  belongs  to  the  E^  component  of  the  mutually  scattered  field 
from  cylinder  and  reflects  for  the  vector  formulation  developed  here. 
Therefore 


-jkb  cos(0  -<j>  )  y  J  <ka)  3 

2Es  -  -jE  •  *  L  -JL -  „(2)  (kR  )  e 

°  n  H  2  (ka)  n 

n 


„  J  (ka)  -j®$. 

£  <-j)  - e  (2kb)  e 

£  H^2  (ka) 

m 


■j  (n-m)0 


2HS  »  H  e 
-  z  o 


■jkb  cos(0_+-4>J  ^ 


r  JA(ka)  ii\  3n<t>. 

L  -S- -  H  2  (kR  )  e 

n  H*  2  (ka)  n 
n 


J’  (ka)  -jnvh . 


-j(n-m)0  , 


5V# 


?■:<$ 

■r. >. 

yv\->, 

i*  4  A* : 

r>>;. 

I  ^ 

W 

vt\W 

II 


nV-*. 


As  kR»l,  the  range-corrected  scattered  far  field  from  the  cylinder 
due  to  the  incoming  wave  from  the  "+"  cylinder  for  each  linear  polariza¬ 
tion  will  be. 
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By  interchanging  the  subscripts  "+"  and  in  eqs.  (38)  and  (39),  we  can 

2  s  2  s 

obtain  the  mutually  scattered  far  fields  +E^  and  from  "+"  cylinder 
due  to  the  incoming  wave  or  directly  scattered  from  the  cylinder. 

The  mutually  scattered  field  from  two  cylinders  for  each  linear  polariza¬ 
tion  is  the  summation  of  these  two  terms  i.e.. 


.Vv  -Vs-VV \ 
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From  which  the  circular  polarization  expression, 

V  «  (V  +  j2ES)//2  (42) 

R  <j>  J  Z 

and 

2E?  -  (2E "  -j  2ES)//2  (43) 

L  <p  z 

can  be  obtained.  Also  here,  because  of  the  unequal  magnitude  and  phase  angle 
of  the  RHCP  and  LHCP  scattered  components  in  the  above  equations,  the  mutually 

scattered  far  field  2E8  from  two  conducting  cylinders  illuminated  with 
RHCP  plane  wave  is  elliptically  polarized. 


It  is  interesting  to  note  that  egs.  (38)  and  (39)  have  leading  phase 


terms  exp ( j2kbsin$ .  sin8  , ) .  This  implies  that  contributions  to  the  image 
2  s  2  s 

reconstructed  from  _E  (or  +E  )  will  have  an  additional  displacement  bsin^K 
(or  -bsin<fc  J  in  the  y-direction.  The  resulting  image  detail  can  be  useful 

therefore  in  some  applications  for  determining  the  bistatic  scattering  angle , 
if  it  is  not  known,  from  the  received  scattered  far  field  of  a  two-cylinder  object 
It  can  be  seen  that  the  second  order  scattering  coefficients  consist 
of  sums  of  products  of  two  scattering  coefficients  of  the  single  cylinder 
and  terms  depending  on  the  distance  between  the  two  cylinders.  This  mutually 
scattered  far  field  expression  can  be  then  physically  interpreted  as 
arising  from  two  scattering  processes  as  the  incident  plane  wave  is  doubly 
scattered  by  one  cylinder  onto  the  other  and  back.  The  second  order  scat¬ 
tering  coefficient  will  approach  zero  for  large  separation  between  the  two 

(2) 

cylinders  due  to  the  behavior  of  the  Hankel  function  H  (2kb) .  This  im- 

n— in 

plies  the  higher  order  scattering  coefficients  could  be  neglected  when  the 
cylinder  spacing  is  large  compared  to  the  illuminating  wavelength. 

In  the  next  section  we  will  examine  the  mutual  scattering  of  two 
cylinders  for  each  polarization  using  numerical  simulation. 


D .  Numerical  Results 

2  s 

Shown  in  Pig.  7  is  the  mutually  scattered  field  components  E^  and 

2E3  calculated  from  eqs.  (42)  and  (43)  for  aspect  angle  9  *  9Cf  with  n  and 

R  _+ 

m  in  the  sunanations  ranging  from  -16  to  16.  The  magnitude  of  the  depolarized 
component  E3  is  seen  to  be  up  to  about  15dB  greater  than  that  of  E^  due  to 
the  effect  of  the  mutual  scattering  between  the  two  cylinders.  The  corres¬ 
ponding  range  profiles  reconstructed  from  each  polarization  component  are 


given  in  Fig.  7(c).  Each  of  these  shows  two  primary  peaks  whose  separation 


is  related  to  the  bistatic  angle  . .  In  addition  the  profile  derived  from 
2  s 

the  El  component  shows  additional  weaker  peaks  and  also  the  higher  order  terms 

Xj 

ignored  in  our  analysis.  The  value  of  i^.  =  6°  is  used  and  those  of  a  and  b 
are  as  before. 

The  numerically  generated  Fourier  space  slices  obtained  by  computing 
2  s  2  s 

E  and  E  as  a  function  of  frequency  and  rotation  angle  0  are  shown  in 
R  L  *t 

Fig.  8.  The  same  polar  format  used  earlier  i.e.,  64  equally  spaced  fre¬ 
quency  points  covering  the  range  (6-17)  GHz  per  radial  line  and  256  radial  lines 
covering  an  angular  range  of  360°,  is  utilized.  The  images  obtained  by 
Fourier  inversion  of  these  slices  are  shown  in  Fig.  9.  The  reconstructed 
images  exhibit  the  effect  of  mutual  scattering  between  two  cylinders. 

The  effect  consists  of  two  symmetrical  "smile  patterns"  with  width  2 (b-a) . 

The  two  peaks  presented  in  Fig.  7  correspond  to  a  cross-sectional  view  of 
the  "smile  pattern"  at  0  ■  90°.  The  lip  separation,  given  by  bsin<f> . , 

is  determined  by  the  bistatic  angle.  It  becomes  wider  as  the  bistatic 
angle  increases  as  indicated  in  Fig.  9(c)  for  which  the  bistatic  angle 
<1^  is  8*  as  compared  to  <|k  »  6°  in  (a)  and  (b)  of  the  figure.  The  smile 
pattern  will  merge  into  a  "ring  pattern"  as  shown  in  Fig.  9(d)  when 

*  0°  i.e.,  the  case  of  monostatic  scattering.  This  suggests  that  the 
bistatic  angle  spanned  by  transmitting  and  receiving  antennas  can,  if 
necessary,  in  certain  situations  be  determined  from  the  lip  separation  of 
the  cross-polarized  image  of  a  two  cylinder  test  target. 

The  total  scattered  wavefield  from  two  cylinders  is  the  superposition 
of  directly  scattered  components  given  in  eqs.  (23),  (24)  and  the  mutually 
scattered  components  given  in  eqs.  (42),  (43)  for  each  state  of  polarization. 
Fourier  space  slices,  obtained  from  the  computed  total  RHCP  and  LHCP 


components,  and  their  reconstructed  images  are  shown  in  Fig.  10.  Due  to 

the  pronounced  depolarization  effect  contributed  from  the  mutual  scattering, 

2  s 

the  RHCP  component  dominated  by  the  E  wavefield,  is  about  up  to  15dB 

K 

1  S 

greater  them  E^.  Therefore  the  reconstructed  image  consists  of  two 

distinct  image  of  the  smile  pattern  bridging  two  adjacent  faint  circles 

representing  the  projection  image  of  the  two  cylinders.  The  slight  spacing 

at  the  corners  of  the  "smile  pattern",  as  compared  to  the  narrower  spacing 

in  the  experimental  results  presented  in  the  following  section,  is  attributed 

to  the  finite  number  of  terms  taken  in  the  calculation  of  eqs.  (40)  and  (41) . 

The  LHCP  component  gives  the  same  image  as  Fig.  6(a)  of  the  directly  scattered 

data  which  is  merely  the  projection  of  the  two  cylinders  because  for  this  case 
2  s  X  s 

E  is  more  them  20dB  weaker  than  E„  and  hence  this  image  is  seen  to  be  free 
L  L 

of  any  depolarization  effects.  Note  in  the  strict  sense  of  imaging,  the  "smile 
patterns"  in  Fig.  9  do  not  correspond  to  any  physical  detail  of  the  two  cylinders 
and  hence  represent  "ghost  images"  or  image  artefact. 

IV.  EXPERIMENTAL  RESULTS 

Due  to  continuing  advances  in  microwave  instrumentation  and  technology, 
it  is  now  possible  to  develop  efficient  and  highly  accurate  computer- 
controlled  microwave  measurement  system  capable  of  collecting  full  vector 
information  of  the  scattered  wavefield  including  amplitude,  phase,  fre¬ 
quency  and  polarization  [1-4],  The  experimental  arrangement  used  to  ac¬ 
quire  wideband  polarization-selective  scattering  data  of  test  objects  is 
shown  in  Fig.  11.  The  test  object,  consisting  of  two  metallic  cylinders 
27.4  cm  apart,  each  100  cm  long  and  3.8  cm  in  radius,  is  mounted  on  a 
computer  controlled  positioner  situated  in  an  anechoic  chamber.  The 


entire  measurement  sequence  consisting  of  angular  positioning  of  the  test 
object  in  azimuth  in  increments  of  1.4®,  incremental  frequency  stepping 
over  a  (6-17)  GHz  spectral  range,  switching  of  antenna  polarization, 
digitization  and  storage  of  the  amplitude  and  phase  of  the  scattered  field 
measured  by  the  coherent  receiver, is  carried  out  under  computer  control 
which  also  performs  the  required  digital  signal  processing.  As  mentioned 
earlier  in  Section  II,  a  single  metallic  cylinder  is  used  to  characterize  the 

system  frequency  response  including  determination  of  the  range  information 
needed  to  correct  for  the  range  term.  Further  detail  on  this  and  on  the  re¬ 
moval  of  clutter,  signals  in  anechoic  chamber  measurements  are  found  in  [3,4]. 

Results  of  two  circular  polarization  measurements,  shown  in  Fig.  12, 
are  obtained  with  the  system  using  RHCP  plane  wave  illumination  produced  by  a 
parabolic  reflector  antenna  and  an  electronically  switchable  RHCP  and  I.HCP  dual¬ 
polarization  horn  receiving  antenna.  The  Fourier  slice  data  and  corres¬ 
ponding  reconstructed  images  for  each  measured  polarization  are  shown. 

Several  observations  can  be  made.  First,  both  the  RHCP  and  LHCP  scattered 
fields  and  their  reconstructed  images  agree  well  with  the  analytical 
results  given  in  Fig.  10.  Secondly,  the  RHCP  component  gives  a  dis¬ 
tinctive  image  of  the  "smile  pattern"  and  two  faint  circles.  This  experi¬ 
mental  result  shows  that  the  higher  order  scattering  terms,  beyond  -16<n<16 
used  in  our  computations,  are  negligibly  small  when  the  cylinders  spacing 

is  chosen  to  be  greater  than  the  radius  of  each  cylinder  and  than  the  illuminating 

2  s 

wavelength.  Thirdly,  the  depolarized  image  of  is  nearly  not  visible 

la 

because  it  is  40db  lower  in  image  intensity  (corresponding  to  20dB  weaker 

1  s 

amplitude)  than  the  F  component . 
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V. 


DISCUSSION 


In  an  attempt  to  clarify  the  role  of  polarization  diversity  in  micro- 
wave  image  formation,  a  systematic  analytical  and  experimental  study  of 
the  directly  and  mutually  scattered  wavefields  of  two  cylinders  in  an 
orthogonal  polarization  basis  was  undertaken .  The  results  show  that  the 
analytical  and  experimental  results  obtained  and  reported  here  are  in  good 
agreement.  These  results  suggest  that  in  the  case  of  dual-cylinder  objects 
illuminated  with  a  RHCP  plane  wave,  the  LHCP  scattered  field  primarily 
determined  by  specular  reflections  from  the  illuminated  surfaces  of 
cylinders  is  free  of  any  noticeable  depolarization  effects  arising  from 
multiple  scattering  between  the  two  cylinders.  Two-dimensional  Fourier 
inversion  of  the  LHCP  Fourier  slice  data  gives  a  projection  image  showing 
the  shape  and  location  of  scattering  centers  on  the  two  cylinders.  Thus 
RHCP  plane  wave  illumination  and  LHCP  scattered  field  measurement  can  be 
used  to  suppress  multiple  scattering  effects,  and  therefore  artefacts,  in 
microwave  imaging  of  cylinders  and  other  more  complex  multiple  segment  objects 
The  RHCP  measurement,  on  the  other  hand,  gives  the  mutual  scattering  informa¬ 
tion.  The  image  information  in  this  case  is  an  "artefact"  as  it  does  not 
correspond  to  any  actual  physical  detail  or  optical  image  of  the  object.  The 
mutual  scattering  information  carries  further  detail  about  the  scatterers  and 
the  scattering  geometry.  For  example  as  shown  here  the  artefact  image  detail 
can  be  related  to  the  bistatic  angle  and  the  cylinder  spacing. 

The  results  presented  here  lead  to  better  understanding  of  the  role 
of  polarization  in  microwave  imaging  in  conveying  more  information  about 
multiple  simple-body  scatterers  e.g.,  cylinders  considered  here,  and  more  com¬ 
plex  objects  in  general.  The  methods  presented  also  serve  to  illustrate 
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ABSTRACT 


Microwave  diversity  imaging  techniques  in  which  spectral  angular  (aspect 
dependent)  and  polarization  degrees  of  freedom  are  employed  in  image  forma¬ 
tion  are  shown  to  be  effective  in  high  resolution  projection  imaging  of 
targets  in  the  presence  of  severe  clutter.  Principles,  methodologies,  and 
experimental  results  demonstrating  this  capability  are  presented.  Implica¬ 
tions  of  the  results  for  air-borne  or  space-borne  down  looking  Imaging 
radars  and  imaging  of  internal  structure  of  penetrable  objects  are  also 


discussed. 


I.  INTRODUCTION 


Frequently  microwave  Imaging  situations  arise  where  the  object  being 
imaged  is  situated  against  a  strongly  reflecting  background  that  produces 
severe  clutter  or  is  obstructed  by  clutter  producing  surfaces  or  foreground. 
Examples  are  found  in  air-borne  or  space-borne  down-looking  imaging  radars 
and  systems  for  the  imaging  of  internal  structure  of  penetrable  bodies.  The 
aim  of  this  paper  is  to  demonstrate  the  capability  of  a  high  resolution 
microwave  diversity  imaging  system  in  discriminating  objects  of  interest  in 
the  presence  of  such  severe  clutter.  The  approach  is  based  on  a  bpot-Ligkt 
imaging  mode  in  which  spectral,  angular,  and  polarization  diversities  are 
utilized  to  access  the  three-dimensional  Fourier  space  of  a  coherently 
illuminated  microwave  scene.  A  brief  theoretical  review  of  the  principle 
is  presented  followed  by  a  description  of  the  microwave  imaging  configura¬ 
tion  employed  in  collecting  realistic  data  of  a  representative  situation: 
that  of  a  low  flying  aircraft.  The  antennas  in  the  imaging  configurations 
studied  are  positioned  as  if  they  were  space-borne  looking  downward  towards 
the  target  and  the  ground.  Only  one  two-dimensional  slice  of  the  three- 
dimensional  Fourier  space  is  obtained;  hence  the  reconstructed  image  is  a 
projective  side-view.  The  reconstructed  images  show  clearly  the  ability  to 
distinguish  the  object  image  from  the  background  clutter  simulated  in  our 
experiments  by  ground  planes  overlayed  by  crumbled  aluminum  foil.  The 
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polarization-enhanced  microwave  Images  obtained  are  comparable  visually  to 
|  high  contrast  side-view  photographs  of  the  scene  taken  with  flood  light 

illumination.  Finally  the  effect  of  relative  motion  and  how  to  account  for 
it  in  the  measurement  is  discussed.  The  results  shown  demonstrate  clearly 
j  the  ability  to  form  recognizable  high  resolution  images  of  complex  scatterers 

in  the  presence  of  severe  clutter. 

II.  THEORETICAL  BACKGROUND 

^  It  has  been  demonstrated  that  by  combining  angular,  spectral,  and  polariza¬ 

tion  diversities  in  a  coherent  imaging  system,  projective  or  tomographic 
images  of  complex  shaped  conducting  bodies  can  be  obtained  with  near  optical 
resolution  [1].  The  principles  on  which  such  a  system  is  based  are  obtained 
from  inverse  scattering  theory  [2,3,4],  which  shows  that  under  the  physical 
optics  and  Bom  approximations,  the  scalar  component  associated  with  one 

A 

polarization  of  the  fields  scattered  in  the  direction  of  the  unit  vector  1R 
and  observed  at  a  distance  R  from  an  object  illuminated  by  a  plane  wave  inci- 

/A 

dent  in  the  direction  of  the  unit  vector  1  is  given  by  [5] 

00 

<MP,R)  -  e~jkS  /  y(r)  ejP’r  dr  (1) 

—  00 

where  y(r)  is  the  object  4catte/u.ng  function  which  is  a  function  of  position¬ 
al  vector  r  referred  to  an  arbitrarily  chosen  reference  origin  in  the  object. 
y(r)  represents  the  3-D  distribution  of  scattering  centers  on  the  target , 


assumed  to  be  nonr-dispersive,  i.e. ,  y  is  independent  of  frequency  cu,  and  as  such 

A  A 

is  closely  related  to  the  surface  contour  of  the  object.  The  vector  p*k(l  -1  ), 

R  i 

/A 

is  the  positional  vector  in  the  Fourier  space.  In  particular,  p*2kl_  in  the  mono 
static  case,  i.e.,  when  transmitter  and  receiver  (T/R)  antennas  are  closely  posi¬ 
tioned.  k«w/c  is  the  wavenumber,  c  being  the  velocity  of  light. 


By  correcting  for  the  range-related  term  in  front  of  the  integral  in  eq. 
(1) ,  we  obtain  a  three-dimensional  Fourier  transform  relationship  between  the 
corrected  data  and  the  object  scattering  function, 

r<?)  -  /  Y<7)ejp’r  d7  -  -j  ^  ejkR  ^(F,R)  (2) 

—00 

This  permits  the  retrieval  of  y(r )  by  Fourier  inversion  of  the  measured  data 
ip  corrected  for  the  range  dependent  factor  preceding  the  integral  in  eq.  (1). 
Techniques  for  determining  R  with  sufficient  accuracy  to  perform  the  range 
correction  and  for  correcting  for  the  frequency  response  of  the  measurement 
system  were  described  elsewhere  [1].  Note  that  the  fT-space  data,  or  F (p) »  can 

/\  A 

be  acquired  over  different  values  of  p  by  varying  1  ,  1  ,  and  k.  AnguZaA  cLL- 

K  1 

A  A 

veAA-ity,  or  angular  degrees  of  freedom  involve  changing  1  and/or  1,  while 

K  1 

varying  k  involves  spectral  degrees  of  freedom  and  is  referred  to  as  Apzc&ULt 
dUveAA-Lty .  The  object  scattering  function  depends  generally  on  polarization, 
which  allows  us  to  exploit  pola/iization  diveAAity  to  further  enrich  the 
accessed  volume  of  the  Fourier  space  with  polarization  related  information 
about  the  object.  Even  though  physical  optics  considerations  predict  no  de¬ 
polarization  in  monostatic  scattering  measurements  some  depolarization  exists 
in  practice.  (Depolarization  effects  become  markedly  noticeable  for  the  bi¬ 
static  case).  In  the  experiments  described  below  the  co-polarized  and  cross- 
polarized  components  of  the  scattered  field  are  measured  for  different  status 
polarization  of  the  incident  illumination  and  the  image  intensities  formed  from 
them  are  added  to  obtain  a  polarization  enhanced  image.  A  more  detailed  dis¬ 
cussion  of  the  role  of  polarization  in  microwave  diversity  imaging  will  be 
given  elsewhere  [5]. 

The  reconstructed  image  quality  depends  on  the  size  of  the  p'-space 
accessed  by  measurement.  The  size  of  p'-space  depends  on  the  values  of  p* 


realized  In  the  measurements  as  Imposed  by  the  spectral  and  angular  windows 
realized  by  the  measurement  gear.  In  order  to  avoid  aliasing  of  the  retrieved 
image  care  must  be  taken  to  choose  the  values  of  p  so  as  to  satisfy  the  Nyquist 
criterion. 

The  Inherently  two-dimensional  nature  of  most  image  display  devices 
leads  to  two  types  of  two-dimensional  presentation- that  are  of  prime  interest. 
One  is  tomogAxiphCc,  or  cross-sectional,  presentation  of  a  three-dimensional 
object.  Given  the  acquired  T(p)  data  in  three-dimensional  p-space,  one  may 
use  the  weighted  pAo j action -itice  the.OA.an  [6],  which  is  the  generalization  of 
the  well-known  pA.ojectCon-4JU.ee.  theoAem  [71,  to  reconstruct  and  display  a 
specific  tomogram  of  the  three-dimensional  object.  The  other  presentation  is 
a  pAjojectioe  image  of  the  three-dimensional  body.  From  the  projection-slice 
theorem,  Che  projective  image  can  be  reconstructed  from  a  central  slice,  or 
meridional  slice  of  p-space  data.  The  imaging  geometry  to  be  described  in  the 
next  section  accesses  one  such  meridonal  slice,  hence  the  retrieved  images 
presented  here  are  projective. 


Figure  1  depicts  the  situation  to  be  simulated,  in  which  a  space-borne 
imaging  platform  containing  a  coherent  broad-band  T/R  system  is  looking  down¬ 
ward  toward  an  aircraft  flying  at  low  altitude.  The  imaging  platform  tracks 
the  aircraft  in  a  ipot-Ught  fashion.  For  simplicity,  we  assume  that  the 
platform  moves  relatively  much  faster  than  the  aircraft,  and  that  during  the 
period  in  which  the  platform  acquires  data  for  one  aspect  look  of  the  air¬ 
craft,  the  change  in  aircraft  is  negligible.  In  actual  situations,  this 
assumption  is  not  valid.  Doppler  shift  caused  by  the  variation  of  distance 
between  the  platform  and  the  aircraft  during  the  data  acquisition  period  of 
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Fig.  1.  Simulated  imaging  geometry. 


one  look  complicates  the  data  acquisition  as  well  as  data  processing.  What 
we  simulate  here  is  an  ideal  case  which  removes  the  unwanted  complexity  to 
show  what  can  ultimately  be  achieved  with  an  imaging  system  that  can  account 
for  the  doppler  shift  and  hence  can  access  the  3-D  Fourier  space  of  target  in 
the  presence  of  relative  motion.  Also  shown  in  Fig.  1  is  the  p-space  data 
accessed,  where  Ap  *  2Ak  *  2Aoj/c ,  is  proportional  to  the  spectral  window 
Aw/2tt  utilized.  Note  that  the  data  lie  on  a  polar  grid.  It  represents  a 
meridional  slice  of  the  three-dimensional  p-space  or  Fourier  space  of  the  ob¬ 
ject.  Thus  the  retrieved  images  will  be  projective,  with  projection  direction 
perpendicular  to  the  plane  of  the  p-space  slice. 

To  simulate  strong  ground  clutter,  a  light  cardboard  plate  covered  with 
a  layer  of  crumbled  aluminum  foil  is  employed  as  a  roughened  ground  plate  or 
background.  Two  such  plates  are  used;  one  with  rms  height  variation  greater 
than  the  mean  wavelength  of  the  spectral  window  employed,  the  other  with  rms 


pression  may  provide  both  the  power  and  spectral  window  width  needed.  But  this 
approach  can  be  restricted  by  the  detector  bandwidth,  which  is  typically  a  frac¬ 
tion  of  the  central  frequency.  The  availability  of  computer -tunable  wide 
band  RF  sources,  e.g.,  sweep  oscillators  and  synthesizers,  makes  it  feasible 
to  generate  effectively  an  extremely  narrow  pulse  using  frequency  synthesis.' 

This  indirect  pulse  synthesis  is  the  approach  we  take  to  achieve  high  resolution. 
Typically  a  spectral  window  of  Af  -  11  GHz  covering  a  frequency  range  from  6.1 
to  17.5  GHz  is  used.  The  number  of  frequency  increments  used  to  cover  this 
spectral  range  was  128.  This  provides  a  range  resolution  6R  *  c/2Af.  of  better 
than  1.5  cm.  The  angular  window  A<J) ,  seen  in  Fig.  1,  used  in  acquiring  data  in 
the  results  presented  below  covered  90  degrees.  The  number  of  angular  samples 
or  looks  in  each  window  is  128. 

Data  for  four  possible  combinations  of  circular  polarization  states  of 
the  T/R  antennas  are  collected.  These  are  designated  in  accordance  to  the 
receiver  and  transmitter  polarizations  following  the  nomenclature  used  in  [8]. 

For  example,  RL  represents  the  data  collected  by  a  right  hand  circularly 
polarized  (RHCP)  receiving  antenna  when  the  object  is  illuminated  by  a  left 
hand  circularly  polarized  (LHCP)  plane  wave  produced  by  a  LHCP  transmitting 
antenna.  When  the  T/R  antennas  are  of  opposite  polarization  states,  we  call 
the  measurement  co-polarlzed,  because  of  the  phase  reversal  of  the  incident 
wave  upon  reflection  from  a  metalized  object.  On  the  other  hand,  and  for  the 
same  reason,  cross-polarization  applies  to  measurement  when  the  T/R  antennas 
have  identical  polarization  states.  Polarization  enhanced  images  are  obtained 
by  Incoherent  combination  of  the  Images  retrieved  from  co-  and  cross-polarized 
Fourier  space  data. 


nterpolated  p-space  data  (left)  and  Image  (right) 
f  the  space  shuttle.  Co-polarized  or  LR  case 
R-LHCP,  T-RHCP). 


Fig.  5.  P -space  data  (left)  and  Image  (right)  of  the  shuttle  above  the 
rougher  ground  plate.  Cross-polarized  case  (R-RHCP,  T-RHCP) . 


The  polarization  enhanced  image,  the  result  of  adding  the  image  inten¬ 
sities  of  Figs.  4  and  5,  is  shown  in  Fig.  6  together  with  a  high  contrast  pic¬ 
torial  side  view  of  the  space  shuttle  on  its  mount  shown  for  comparison  taken 
under  flood -light  conditions  to  emphasize  similar  outline.  The  similarity  of 
the  two  Images  is  striking.  Most  of  the  features  of  the  shuttle  are  present 
without  much  degradation  by  the  influence  of  the  clutter  plate.  This  is  attri¬ 
buted  primarily  to  the  wide  dynamic  range  of  the  measurement  system  estimated 
to  be  68  dB,  which  enables  the  delineation  of  weak  scattering  detail  in  the 
presence  of  strong  scatterers,  and  due  to  the  excellent  range  resolution  (~1.5  cm) 
which  enables  isolating  the  clutter  image  spatially. 


Fig.  6.  Polarization  enhanced  image  (RR+LR)  (left)  and  optical  image 
(right). 

Enhanced  images  of  the  shuttle  with  each  of  the  two  clutter  plates  (of 
different  rms  height  variations)  are  shown  in  Fig.  7.  The  image  with  the 
smoother  clutter  plate  compares  less  favorably  to  that  with  rougher  clutter 
plate.  The  shuttle  is  less  visible  and  artifact  is  present  in  the  clutter 
image.  The  weaker  shuttle  Image  is  believed  to  be  due  to  the  less  diffuse 
nature  of  scattering  from  the  smoother  plate  which  tends  to  produce  stronger 
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We  have  demonstrated  for  the  first  time  that  despite  the  presence  of  the 
relatively  strong  background  clutter,  images  of  high  quality  can  be  obtained 
with  a  high  resolution  system.  The  microwave  images  are  comparable  visually 
to  their  visual  counterparts  under  mono-directional  floodlight  illumination. 
Although  for  the  images  shown  here  a  90  degree  angular  window  is  utilized, 
reduction  of  the  window  size  down  to  30  degrees  continued  to  yield  discerni¬ 
ble,  though  somewhat  degraded,  images  of  the  target. 

Separation  of  the  target  from  the  clutter  may  be  desirable.  Spatial 
filtering  in  the  image  plane  can  be  employed  to  achieve  this.  Coherent  sub¬ 
traction  in  the  p-space  is  another  way  in  which  the  p-space  data  of  the 
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tered  noticeably,  which  makes  the  coherent  removal  of  the  clutter  difficult.  In 
actual  situations,  this  strategy  is  impractical  since  it  requires  measurement 
of  the  same  clutter  at  exactly  the  same  aspect  looks  twice;  once  when  the  target 
is  present,  the  other  without.  The  more  feasible  practice  as  stated  earlier 
is  separating  the  clutter  in  the  image  space  by  spatial  filtering.  The  one- 
dimensional  FFT  of  the  data  at  each  look  is  the  range  profile  of  the  scene 
at  that  aspect  look.  With  good  range  resolution,  detail  belonging  to  the 
image  of  the  object  and  that  of  the  clutter  can  be  separated. 

As  mentioned  earlier,  the  simulated  conditions  in  the  above  study  are 
ideal,  l.e.,  the  target  motion  is  Ignored.  In  actuality,  the  Imaging  platform 
and  the  target  are  in  relative  motion.  It  is  however  the  change  of  aspect 
of  the  object  scene  during  motion  that  provides  image  formation;  i.e.,  it  is 
the  changes  of  their  relative  orientations  that  contributes  to  image  formation. 
This  has  bden  the  central  process  in  the  spot-light  imaging  scheme  studied 
here.  The  effect  of  relative  motion  is  merely  to  introduce  doppler  shift  in  the 
frequency  of  the  received  echo.  Compensating  for  the  doppler  shift  in  the 
measurement  methodology  described  here  would  allow  accessing  the  p-space  of 
a  moving  target  in  the  manner  described  above.  Methods  for  removing  the 

doppler  shift  caused  by  relative  motion  are  currently  under  study.  One  pro¬ 
mising  method  being  studied  utilizes  what  we  call  "orchestrated  frequency 
response  measurement"  in  which  the  doppler  shift  of  the  echo  relative  to  the 


transmitted  signal  is  rapidly  determined  and  the  information  used  to  set  the 
frequency  of  a  synthesized  source  used  as  local  oscillator  for  the  receiver 
to  a  value  needed  to  compensate  for  doppler  shift.  The  L.O.  synthesizer  in 
this  scheme  is  phase-locked  to  a  transmitter  synthesizer  providing  the  object 
illumination <  Equipment  for  performing  this  study  obtained  under  a  $350K 
DoD/URIP  grant  is  presently  being  assembled.  The  result  of  this  study  will 
be  reported  on  in  the  near  future. 


VI.  ACKNOWLEDGEMENT 

The  work  reported  here  was  supported  by  a  grant  from  General  Electric 
Co.  and  under  partial  support  from  the  Air  Force  Office  of  Scientific 

Research. 


VII.  REFERENCES 

1.  N.H.  Farhat,  et.  al.,  "Prospects  for  Three-dimensional  Projective 
and  Tomographic  Imaging  Radar  Networks,"  Radio  Science,  Vol.  19, 

No.  5,  pp.  1347-1355,  Sept. -Oct.  1984.  (see  also  Appendix  I  of  this 
document) 

2.  R.M.  Lewis,  "Physical  Optics  Inverse  Diffraction,"  IEEE  Trans,  on 
Antennas  &  Propagations,  Vol.  AP-17,  pp.  308-314,  1969. 

3.  N.  Bo jar ski,  "Inverse  Scattering,"  Final  Report  N000  19-73-C-0312F, 
Nav.  Air  Syst.  Command,  Warminister,  PA.,  February  1974. 

4.  S.R.  Raz,  "On  Scatterer  Reconstruction  from  Far-field  Data," 

IEEE  Trans,  on  Antennas  &  Propagations,  Vol.  AP-24,  pp.  66-70,  1976. 

5.  T.H.  Chu  and  N.H.  Farhat,  "Polarization  Effects  in  Microwave  Diversity 
Imaging  of  Conducting  Cylinders,"  submitted  for  publication. 

6.  C.K.  Chan  and  N.H.  Farhat,  "Frequency  Swept  Tomographic  Imaging  of 
Three-Dimensional  Perfectly  Conducting  Objects,"  IEEE  Trans,  on 
Antennas  &  Propagations.  Vol.  AP-29,  pp.  321-319,  1981. 

7.  R.M.  Mersereau  and  A.V.  Oppenhelm,  "Digital  Reconstruction  of 
Multidimensional  Signals  from  Their  Projections,"  Proc.  IEEE,  Vol. 

62,  pp.  1319-1338,  1978. 

8.  G.T.  Ruck,  et.  al..  Chapter  2  in  Radar  Cross  Section  Handbook,  Vol. 

1,  G.T.  Ruck  (Ed.),  Plenum  Press,  New  York,  1970. 


APPENDIX  III 


Reprinted  from  Applied  Optica,  Vol.  24,  page  14(1,  Bley  IS,  1*85 
Copyright  ©  1985  by  the  Optical  Society  of  America  and  reprinted  by  permiaaion  of  the  copyright  owner. 

Optical  implementation  of  the  Hopfield  model 

Nabil  H.  Farhat,  Demetri  Psaltis,  Aluizio  Praia,  and  Eung  Paek 


Optical  implementation  of  content  addressable  associative  memory  based  on  the  Hopfield  model  for  neural 
networks  and  on  the  addition  of  nonlinear  iterative  feedback  to  a  vector-matrix  multiplier  is  described.  Nu¬ 
merical  and  experimental  results  presented  show  that  the  approach  is  capable  of  introducing  accuracy  and 
robustness  to  optical  processing  while  maintaining  the  traditional  advantages  of  optics,  namely,  parallelism 
and  massive  interconnection  capability.  Moreover  a  potentially  useful  link  between  neural  processing  and 
optics  that  can  be  of  interest  in  pattern  recognition  and  machine  vision  is  established. 


L  Introduction 

It  is  well  known  that  neural  networks  in  the  eye-brain 
system  process  information  in  parallel  with  the  aid  of 
large  numbers  of  simple  interconnected  processing  el¬ 
ements,  the  neurons.  It  is  also  known  that  the  system 
is  very  adept  at  recognition  and  recall  from  partial  in¬ 
formation  and  has  remarkable  error  correction  capa¬ 
bilities. 

Recently  Hopfield  described  a  simple  model1  for  the 
operation  of  neural  networks.  The  action  of  individual 
neurons  is  modeled  as  a  thresholding  operation  and 
information  is  stored  in  the  interconnections  among  the 
neurons.  Computation  is  performed  by  setting  the 
state  (on  or  off)  of  some  of  the  neurons  according  to  an 
external  stimulus  and,  with  the  interconnections  set 
according  to  the  recipe  that  Hopfield  prescribed,  the 
state  of  all  neurons  that  are  interconnected  to  those  that 
are  externally  stimulated  spontaneously  converges  to 
the  stored  pattern  that  is  most  similar  to  the  external 
input.  The  basic  operation  performed  is  a  nearest- 
neighbor  search,  a  fundamental  operation  for  pattern 
recognition,  associative  memory,  and  error  correction. 
A  remarkable  property  of  the  model  is  that  powerful 
global  computation  is  performed  with  very  simple, 
identical  logic  elements  (the  neurons).  The  intercon¬ 
nections  provide  the  computation  power  to  these  simple 
logic  elements  and  also  enhance  dramatically  the  stor- 
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age  capacity;  approximately  N/4  In N  bits/neuron  can 
be  stored  in  a  network  in  which  each  neuron  is  con¬ 
nected  to  N  others.2  Another  important  feature  is  that 
synchronization  among  the  parallel  computing  elements 
is  not  required,  making  concurrent,  distributed  pro¬ 
cessing  feasible  in  a  massively  parallel  structure.  Fi¬ 
nally,  the  model  is  insensitive  to  local  imperfections 
such  as  variations  in  the  threshold  level  of  individual 
neurons  or  the  weights  of  the  interconnections. 

Given  these  characteristics  we  were  motivated  to 
investigate  the  feasibility  of  implementing  optical  in¬ 
formation  processing  and  storage  systems  that  are  based 
on  this  and  other  similar  models  of  associative  memo¬ 
ry.3,4  Optical  techniques  offer  an  effective  means  for 
the  implementation  of  programmable  global  intercon¬ 
nections  of  very  large  numbers  of  identical  parallel  logic 
elements.  In  addition,  emerging  optical  technologies 
such  as  2-D  spatial  light  modulators,  optical  bistability, 
and  thin-film  optical  amplifiers  appear  to  be  very  well 
suited  for  performing  the  thresholding  operation  that 
is  necessary  for  the  implementation  of  the  model. 

The  principle  of  the  Hopfield  model  and  its  impli¬ 
cations  in  optical  information  processing  have  been 
discussed  earlier.6,6  Here  we  review  briefly  the  main 
features  of  the  model,  give  as  an  example  the  results  of 
a  numerical  simulation,  describe  schemes  for  its  optical 
implementation,  then  present  experimental  results 
obtained  with  one  of  the  schemes  and  discuss  their 
implications  as  a  content  addressable  associative 
memory  (CAM). 

II.  Hopfield  Model 

Given  a  set  of  M  bipolar,  binary  (1,-1)  vectors  vi’n), 
i  m  1,2,3 . .  JV,  m  *  1,2,3  . .  .M,  these  are  stored  in  a 
synaptic  matrix  in  accordance  with  the  recipe 

Ti>  -  I  v|""v}"»,  ij  -  1,2,3. .  .N,  T«  -  0,  (1) 

1 

v'm>  are  referred  to  as  the  nominal  state  vectors  of  the 
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memory.  If  the  memory  is  addressed  by  multiplying 
the  matrix  Ty  with  one  of  the  state  vectors,  say  v1(mo),  it 
yields  the  estimate 


where 


*(—>  «  £  Ty  V}™* 

J 

-  £  L  v!">vj">v}*°> 
jw*i  m 

«(N-  l)v‘”“»+  Z 


t«) 


(2) 


(3) 


❖J'"0'  consists  of  the  sum  of  two  terms:  the  first  is  the 
input  vector  amplified  by  (N  —  1);  the  second  is  a  linear 
combination  of  the  remaining  stored  vectors  and  it 
represents  an  unwanted  cross-talk  term.  The  value  of 
the  coefficients  amjno  is  equal  to  y/N  -  1  on  the  average 
(the  standard  deviation  of  the  sum  of  N  -  1  random 
bits),  and  since  (M  —  1)  such  coefficients  are  randomly 
added,  the  value  of  the  second  term  will  on  the  average 
be  equal  to  V(M  -  l)(N  -  1).  If  N  is  sufficiently 
larger  than  M,  with  high  probability  the  elements  of  the 
vector  ♦}mo)  will  be  positive  if  the  corresponding  ele¬ 
ments  of  v‘mo>  are  equal  to  +1  and  negative  otherwise. 
Thresholding  of  ❖}n,ffl)  will  therefore  yield  v}moi: 


vj"">  -  sgn[*‘,“)]  ■ 


+1  if  >  0 
— 1  otherwise. 


(4) 


When  the  memory  is  addressed  with  a  binary  valued 
vector  that  is  not  one  of  the  stored  words,  the  vector- 
matrix  multiplication  and  thresholding  operation  yield 
an  output  binary  valued  vector  which,  in  general,  is  an 
approximation  of  the  stored  word  that  is  at  the  shortest 
Hamming  distance  from  the  input  vector.  If  this  out¬ 
put  vector  is  fed  back  and  used  as  the  input  to  the 
memory,  the  new  output  is  generally  a  more  accurate 
version  of  the  stored  word  and  continued  iteration 
converges  to  the  correct  vector. 

The  insertion  and  readout  of  memories  described 
above  are  depicted  schematically  in  Fig.  1.  Note  that 
in  Fig.  1(b)  the  estimate  $ jma>  can  be  viewed  as  the 
weighted  projection  of  Ty.  Recognition  of  an  input 
vector  that  corresponds  to  one  of  the  state  vectors  of  the 
memory  or  is  close  to  it  (in  the  Hamming  sense)  is 
manifested  by  a  stable  state  of  the  system.  In  practice 
unipolar  binary  (0,1)  vectors  or  words  b{m>  of  bit  length 
N  may  be  of  interest.  The  above  equations  are  then 
applicable  with  [2b,!B,)  -  1)  replacing  vjm>  in  Eq.  (1)  and 
bfmo>  replacing  v,tmo,  in  Eq.  (2).  For  such  vectors  the 
SNR  of  the  estimate  ♦}mo)  can  be  shown  to  be  lower  by 
a  factor  of  y/2.1 

An  example  of  the  Ty  matrix  formed  from  four  binary 
unipolar  vectors,  each  being  IV  *  20  bits  long,  is  given 
in  Fig.  2  along  with  the  result  of  a  numerical  simulation 
of  the  process  of  initializing  the  memory  matrix  with  a 
partial  version  of  b}4)  in  which  the  first  eight  digits  of  b}4> 
are  retained  and  the  remainder  set  to  zero.  The 
Hamming  distance  between  the  initializing  vector  and 
b}4>  is  6  bits  and  it  is  9  or  more  bits  for  the  other  three 
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(a) 
Fig.  1. 


READ-IN  (b)  READ-OUT 

(a)  Insertion  and  (b)  readout  of  memories. 


stored  vectors.  It  is  seen  that  the  partial  input  is  rec¬ 
ognized  as  b]4)  in  the  third  iteration  and  the  output  re¬ 
mains  stable  as  b|4)  thereafter.  This  convergence  to  a 
stable  state  generally  persists  even  when  the  Ty  matrix 
is  binarized  or  clipped  by  replacing  negative  elements 
by  minus  ones  and  positive  elements  by  plus  ones  evi¬ 
dencing  the  robustness  of  the  CAM.  A  binary  synaptic 
matrix  has  the  practical  advantage  of  being  more  readily 
implementable  with  fast  programmable  spatial  light 
modulators  (SLM)  with  storage  capability  such  as  the 
Litton  Lightmod.7  Such  a  binary  matrix,  implemented 
photographically,  is  utilized  in  the  optical  implemen¬ 
tation  described  in  Sec.  Ill  and  evaluated  in  Sec.  IV  of 
this  paper. 

Several  schemes  for  optical  implementation  of  a  CAM 
based  on  the  Hopfield  model  have  been  described  ear¬ 
lier.5  In  one  of  the  implementations  an  array  of  light 
emitting  diodes  (LEDs)  is  used  to  represent  the  logic 
elements  or  neurons  of  the  network.  Their  state  (on  or 
off)  can  represent  unipolar  binary  vectors  such  as  the 
state  vectors  b,(m)  that  are  stored  in  the  memory  matrix 
Ty.  Global  interconnection  of  the  elements  is  realized 
as  shown  in  Fig.  3(a)  through  the  addition  of  nonlinear 
feedback  (thresholding,  gain,  and  feedback)  to  a  con¬ 
ventional  optical  vector-matrix  multiplier8  in  which  the 
array  of  LEDs  represents  the  input  vector  and  an  array 
of  photodiodes  (PDs)  is  used  to  detect  the  output  vector. 
The  output  is  thresholded  and  fed  back  in  parallel  to 
drive  the  corresponding  elements  of  the  LED  array. 
Multiplication  of  the  input  vector  by  the  Ty  matrix  is 
achieved  by  horizontal  imaging  and  vertical  smearing 
of  the  input  vector  that  is  displayed  by  the  LEDs  on  the 
plane  of  the  Ty  mask  [by  means  of  an  anamorphic  lens 
system  omitted  from  Fig.  3(a)  for  simplicity].  A  second 
anamorphic  lens  system  (also  not  shown)  is  used  to 
collect  the  light  emerging  from  each  row  of  the  Ty  mask 
on  individual  photosites  of  the  PD  array.  A  bipolar  Ty 
matrix  is  realized  in  incoherent  light  by  dividing  each 
row  of  the  Ty  matrix  into  two  subrows,  one  for  positive 
and  one  for  negative  values  and  bringing  the  light 
emerging  from  each  subrow  to  focus  on  two  adjacent 
photosites  of  the  PD  array  that  are  electrically  con¬ 
nected  in  opposition  as  depicted  in  Fig.  3(b).  In  the 
system  shown  in  Fig.  3(a),  feedback  is  achieved  by 
electronic  wiring.  It  is  possible  and  preferable  to  dis¬ 
pose  of  electronic  wiring  altogether  and  replace  it  by 
optical  feedback.  This  can  be  achieved  by  combining 
the  PD  and  LED  arrays  in  a  single  compact  hybrid  or 
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Fig.  3.  Concept  for  optical  implementation  of  a  content  addressable 
memory  baaed  on  the  Hopfield  model,  (a)  Matrix-vector  multiplier 
incorporating  nonlinear  electronic  feedback,  (b)  Scheme  for  realizing 
a  binary  bipolar  memory  mask  transmittance  in  incoherent  light. 


monolithic  structure  that  can  also  be  made  to  contain 
all  ICs  for  thresholding,  amplification,  and  driving  of 
LEDs.  Optical  feedback  becomes  even  more  attractive 
when  we  consider  that  arrays  of  nonlinear  optical  light 
amplifiers  with  internal  feedback9  or  optical  bistability 


devices  (OBDs)10  can  be  used  to  replace  the  PD/LED 
arrays.  This  can  lead  to  simple  compact  CAM  struc¬ 
tures  that  may  be  interconnected  to  perform  higher- 
order  computations  than  the  nearest-neighbor  search 
performed  by  a  single  CAM. 

We  have  assembled  a  simple  optical  system  that  is  a 
variation  of  the  scheme  presented  in  Pig.  3(a)  to  simu¬ 
late  a  network  of  N  =  32  neurons.  The  system,  details 
of  which  are  given  in  Figs.  5-8,  was  constructed  with  an 
array  of  thirty-two  LEDs  and  two  multichannel  silicon 
PD  arrays,  each  consisting  of  thirty-two  elements. 
Twice  as  many  PD  elements  as  LEDs  are  needed  in 
order  to  implement  a  bipolar  memory  mask  transmit¬ 
tance  in  incoherent  light  in  accordance  with  the  scheme 
of  Fig.  3(b).  A  bipolar  binary  T„  mask  was  prepared 
for  M  =  3  binary  state  vectors.  The  three  vectors  or 
words  chosen,  their  Hamming  distances  from  each 
other,  and  the  resulting  T,;  memory  matrix  are  shown 
in  Fig.  4.  The  mean  Hamming  distance  between  the 
three  vectors  is  16.  A  binary  photographic  transpar¬ 
ency  of  32  X  64  square  pixels  was  computer  generated 
from  the  Ty  matrix  by  assigning  the  positive  values  in 
any  given  row  of  Ty  to  transparent  pixels  in  one  subrow 
of  the  mask  and  the  negative  values  to  transparent 
pixels  in  the  adjacent  subrow.  To  insure  that  the  image 
of  the  input  LED  array  is  uniformly  smeared  over  the 
memory  mask  it  was  found  convenient  to  split  the  mask 
in  two  halves,  as  shown  in  Fig.  5,  and  to  use  the  resulting 
submasks  in  two  identical  optical  arms  as  shown  in  Fig. 
6.  The  size  of  the  subrows  of  the  memory  submasks  was 
made  exactly  equal  to  the  element  size  of  the  PD  arrays 
in  the  vertical  direction  which  were  placed  in  register 
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Clipped  aaaory  matrix: 

II  -1  1  1-1-1  1  1  1  1-1  l  1  1  1-1  1-1  1  1  1  1-1-1  l-l  1  -1  -1  -1  -1  -1 

-1  0  1  -1  -1  -1  -1  -1  -1  -1  -1  -1  -1  1-1  1  1  1-1  1  1  1  t  1  1  1  -1  -11-11  -1 

1  1  »  -1  -1  -1  -t  1-1  1  1  1  1  1-1  1-1-1  1-1  1  1-1  1  -1  -1  -1  1-1-1  1-1 

1-1-1  #  1  1  1  1  1  1  1  1  1-1  1  -1  -1  -1  1  -l  -1  -1  -1  -l  -1  -1  1  l-l  1-1  1 

-1-1-1  1  9  1  l-l  1  -1  -1  -1  -1  -1  1-1  1  1-1  1-1-1  1-1  1  1  1-1  1  l-l  1 

-1  -l  -1  1  1  0  1-1  t  -1  -1  -1  -1  -1  1-1  1  1-1  1-1-1  1-1  1  1  1-1  1  1-1  1 

1-1-1  1  1  1  9  1  1  1  1  1  1-1  1  -1  -1  -1  1  -1  -1  -1  -1  -l  -1  -1  1  1-1  1-1  1 

1-1  1  1-1-1  l  0  1  1-1  1  l  1  1-1  1-1  1  1  1  1-1-1  1-1  1  -1  -t  -1  -1  -1 

1-1-1  1  1  1  l  1  9  1  1  1  1-1  1  -1  -1  -1  1  -1  -1  -1  -1  -1  -1  -1  1  1-1  1-1  1 

1-1  1  1-1-1  1  1  l  9-1  1  1  1  l-l  1-1  1  1  1  1-1-1  1-1  1  -1  -t  -1  -1  -1 

-1-1  1  1-1-1  1-1  1-19-1-1  1  1-1-1  1-1-1  1  1  1  -t  -1  1  1  l  1-1-1  -1 
1-1  1  1-1-1  1  1  1  1-1  0  1  1  1-1  1-1  1  1  1  1-1-1  1-1  1  -1  -l  -1  -1  -1 

1-1  1  1-1-1  1  1  1  1-1  1  9  l  1-1  l-l  1  1  1  1-1-1  1-1  1  -l  -1  -1  -l  -1 

1  11-1-1-1-1  1-1  1  1  1  1  9-1  1  -l  -1  1-1  1  1-1  1  -l  -1  -1  1  -1  -1  1  -1 
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1  1  1  -1  -1  -l  -1  1-1  1  1  1  1  1-1  l  -1  -1  1-1  9  1-1  l  -1  -1  -1  1-1-1  1-1 

1  l  1  -1  -1  -1  -1  1-1  1  l  1  1  1-1  1  -l  -l  1-1  l  9-1  1  -1  -l  -1  1  -l  -1  l-l 
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-l  1  -1  -1  1  1  -1  -1  -1  -1  1  -1  -1  -1  -1  1-1  1  -1  -1  -1  -1  1  l-l  9-1  t  1  1  1  1 

l  -1  -1  1  1  1  1  1  1  1  1  1  1-1  1  -l  -1  -1  1  -l  -1  -1  -1  -1  -l  -1  9  1-1  1-1  1 

-1-1  1  1-1-1  1-1  1-1  l  -1  -1  1  1-1-1  1-1-1  l  1  l  -l  -l  1  1  9  1-1  -1  -1 

-1  l  -1  -1  1  1  -1  -1  -1  -1  1  -1  -1  -1  -1  1-1  1  -1  -1  -1  -l  1  1-1  1-1  1  9  1  1  1 

-1-1-1  1  1  1  l-l  l  -l  -1  -1  -1  -1  l-l  l  1-1  l  -1  -1  l-l  1  1  l-l  1  9-1  1  Fig-  *■  Stored  words,  their 
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-l  -l  -l  1  1  1  t  -1  l  -l  -1  -1  -1  -l  1-1  1  l-l  1-1-1  l-l  1  1  l-l  1  1-1  9  clipped  T,;  memory  matrix. 


Fig.  5.  Two  halves  of  T,;  memory  mask. 


against  the  masks.  Light  emerging  from  each  subrow 
of  a  memory  submask  was  collected  (spatially  inte¬ 
grated)  by  one  of  the  vertically  oriented  elements  of  the 
multichannel  PD  array.  In  this  fashion  the  anamorphic 
optics  required  in  the  output  part  of  Fig.  3(a)  are  dis¬ 
posed  of,  resulting  in  a  more  simple  and  compact  sys¬ 
tem.  Pictorial  views  of  the  input  LED  array  and  the 


two  submask/PD  array  assemblies  are  shown  in  Figs. 
7(a)  and  (b),  respectively.  In  Fig.  7(b)  the  left  memory 
submask/PD  array  assembly  is  shown  with  the  submask 
removed  to  reveal  the  silicon  PD  array  situated  behind 
it.  All  electronic  circuits  (amplifiers,  thresholding 
comparators,  LED  drivers,  etc.)  in  the  thirty-two  par¬ 
allel  feedback  channels  are  contained  in  the  electronic 
amplification  and  thresholding  box  shown  in  Fig.  6(a) 
and  in  the  boxes  on  which  the  LED  array  and  the  two 
submask/PD  array  assemblies  are  mounted  (see  Fig.  7). 
A  pictorial  view  of  a  composing  and  display  box  is  shown 
in  Fig.  8.  This  contains  an  arrangement  of  thirty-two 
switches  and  a  thirty-two  element  LED  display  panel 
whose  elements  are  connected  in  parallel  to  the  input 
LED  array.  The  function  of  this  box  is  to  compose  and 
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Fig.  7.  Views  of  (a)  input  LED  array  and  (b)  memory  submask/PD 
array  assemblies. 


display  the  binary  input  word  or  vector  that  appears  on 
the  input  LED  array  of  the  system  shown  in  Fig.  7(a). 
Once  an  input  vector  is  selected  it  appears  displayed  on 
the  composing  box  and  on  the  input  LED  box  simulta¬ 
neously.  A  single  switch  is  then  thrown  to  release  the 
system  into  operation  with  the  composed  vector  as  the 


initializing  vector.  The  final  state  of  the  system,  the 
output,  appears  after  a  few  iterations  displayed  on  the 
input  LED  array  and  the  display  box  simultaneously. 
The  above  procedure  provides  for  convenient  exercising 
of  the  system  in  order  to  study  its  response  vs  stimulus 
behavior.  An  input  vector  is  composed  and  its 
Hamming  distance  from  each  of  the  nominal  state 
vectors  stored  in  the  memory  is  noted.  The  vector  is 
then  used  to  initialize  the  CAM  as  described  above  and 
the  output  vector  representing  the  final  state  of  the 
CAM  appearing,  almost  immediately,  on  the  display  box 
is  noted.  The  response  time  of  the  electronic  feedback 
channels  as  determined  by  the  3-dB  roll-off  of  the  am¬ 
plifiers  was  ~60  msec.  Speed  of  operation  was  not  an 
issue  in  this  study,  and  thus  low  response  time  was 
chosen  to  facilitate  the  experiment. 

IV.  Results 

The  results  of  exercising  and  evaluating  the  perfor¬ 
mance  of  the  system  we  described  in  the  preceding 
section  are  tabulated  in  Table  I.  The  first  run  of  ini¬ 
tializing  vectors  used  in  exercising  the  system  were  error 
laden  versions  of  the  first  word  b|l).  These  were  ob¬ 
tained  from  b}1*  by  successively  altering  (switching)  the 
states  of  1,2,3 ...  up  to  AT  of  its  digits  starting  from  the 
Nth  digit  In  doing  so  the  Hamming  distance  between 
the  initializing  vector  and  bjn  is  increased  linearly  in 
unit  steps  as  shown  in  the  first  column  of  Table  I 
whereas,  on  the  average,  the  Hamming  distance  be¬ 
tween  all  these  initializing  vectors  and  the  other  two 
state  vectors  remained  approximately  the  same,  about 
N/2  =  16.  The  final  states  of  the  memory,  i.e.,  the 
steady-state  vectors  displayed  at  the  output  of  the 
system  (the  composing  and  display  box)  when  the 
memory  is  prompted  by  the  initializing  vectors,  are 
listed  in  column  2  of  Table  I.  When  the  Hamming 
distance  of  the  initializing  vector  from  b)u  is  <11,  the 
input  is  always  recognized  correctly  as  bj1*.  The  CAM 
is  able  therefore  to  recognize  the  input  vector  as  bj1*  even 
when  up  to  11  of  its  digits  (37.5%)  are  wrong.  This 
performance  is  identical  to  the  results  obtained  with  a 
digital  simulation  shown  in  parenthesis  in  column  2  for 
comparison.  When  the  Hamming  distance  is  increased 
further  to  values  lying  between  12  and  22,  the  CAM  is 
confused  and  identifies  erroneously  other  state  vectors, 
mostly  b[3),  as  the  input.  In  this  range,  the  Hamming 
distance  of  the  initializing  vectors  from  any  of  the  stored 
vectors  is  approximately  equal  making  it  more  difficult 
for  the  CAM  to  decide.  Note  that  the  performance  of 
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Table  L  Optical  CAM  Partaonaoca 


Hamming 
distance  of 
initializing 
vector  from 
b‘-> 

Recognized 
vector 
(m  ■  1) 

Recognized 

vector 

(m  »  2) 

Recognized 
vector 
(m  -  3) 

0 

1(1) 

2(2) 

3(3) 

1 

1(1) 

2(2) 

3(3) 

2 

1(1) 

2(2) 

3(3) 

3 

1(1) 

2(2) 

3(3) 

4 

1(1) 

2(2) 

3(3) 

5 

1(1) 

2(2) 

3(3) 

6 

1(1) 

2(2) 

3(3) 

7 

1(1) 

2(2) 

3(3) 

8 

1(1) 

2(2) 

3(3) 

9 

1(1) 

2(2) 

3(3) 

10 

1(1) 

1(1) 

3(3) 

11 

1(1) 

2(2) 

3(3) 

12 

3(3) 

3,2  (3) 

3(3) 

13 

3(3) 

3(3) 

3(3) 

14 

3(3) 

1.3(1) 

3(3) 

15 

1  (OSC) 

l(l) 

2,3  (3) 

16 

3  (OSC) 

1(1) 

3(3) 

17 

3  (OSC) 

1 (OSC) 

3(3) 

18 

3(3) 

1(3) 

3  (OSC) 

19 

3(3) 

2(3) 

3(3) 

20 

3(1) 

3(3) 

3  (OSC) 

21 

1.2(1) 

3(3) 

3  (OSC) 

22 

3<T) 

3(3) 

3  (OSC) 

23 

T(I) 

3(3) 

3  (OSC) 

24 

I  (T) 

3(3) 

3(3) 

25 

T(T) 

3(5) 

3(3) 

26 

Id) 

3(3) 

3(3) 

27 

KT) 

3(3) 

3(3) 

28 

T  (T) 

3(3) 

3(3) 

29 

T(T) 

3(3) 

3(3) 

30 

KT) 

3(3) 

3(3) 

31 

Id) 

2(3) 

3(3) 

32 

id) 

3(3) 

3(3) 

the  CAM  and  results  of  digital  simulation  in  this  range 
of  Hamming  distance  are  comparable  except  for  the 
appearance  of  oscillations  (designated  by  OSC)  in  the 
digital  simulation  when  the  outcome  oscillated  between 
several  vectors  that  were  not  the  nominal  state  vectors 
of  the  CAM.  Beyond  a  Hamming  distance  of  22  both 
the  optical  system  and  the  digital  simulationjdentified 
the  initializing  vectors  as  the  complement  b|l)  of  b}M. 
This  is  expected  because  it  can  be  shown  using  Eq.  (1) 
that  the  Ty  matrix  formed  from  a  set  of  vectors  bjjjj1  is 
identical  to  that  formed  by  the  complementary  set  b1,m). 
The  complementary  vector  can  be  viewed  as  a  contrast 
reversed  version  of  the  original  vector  in  which  zeros  and 
ones  are  interchanged.  Recognition  of  a  complemen¬ 
tary  state  vector  by  the  CAM  is  analogous  to  our  rec¬ 
ognizing  a  photographic  image  from  the  negative. 

Similar  results  of  initializing  the  CAM  with  error 
laden  versions  of  bj®  and  bf3>  were  also  obtained.  These 
are  presented  in  columns  2  and  3  of  Table  I.  Here  again 
we  see  when  the  Hamming  distance  of  the  initializing 
vector  from  b}3),  for  example,  ranged  between  1  and  14, 
the  CAM  recognized  the  input  correctly  as  bf3>  as  shown 
in  column  3  of  the  table  and  as  such  it  did  slightly  better 
than  the  results  of  digital  simulation.  Oscillatory  be¬ 
havior  is  also  observed  here  in  the  digital  simulation 
when  the  range  of  Hamming  distance  between  the  ini- 
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tializing  vector  from  all  stored  vectors  approached  the 
mean  Hamming  distance  between  the  stored  vectors. 
Beyond  this  range  the  memory  recognizes  the  input  as 
the  complementary  of  b|3>. 

In  studying  the  results  presented  in  Table  I  several 
observations  can  be  made:  The  optically  implemented 
CAM  is  working  as  accurately  as  the  digital  simulations 
and  perhaps  better  if  we  consider  the  absence  of  oscil¬ 
lations.  These  are  believed  to  be  suppressed  in  the 
system  because  of  the  nonsharp  thresholding  performed 
by  the  smoothly  varying  nonlinear  transfer  function  of 
electronic  circuits  compared  with  the  sharp  thresh¬ 
olding  in  digital  comptations.  The  smooth  nonlinear 
transfer  function  and  the  finite  time  constant  of  the 
optical  system  provide  a  relaxation  mechanism  that 
substitutes  for  the  role  of  asynchronous  switching  re¬ 
quired  by  the  Hopfield  model.  Generally  the  system 
was  able  to  conduct  successful  nearest-neighbor  search 
when  the  inputs  to  the  system  are  versions  of  the  nom¬ 
inal  state  vectors  containing  up  to  ~30%  error  in  their 
digits.  It  is  worth  noting  that  this  performance  is 
achieved  in  a  system  built  from  off-the-shelf  electronic 
and  optical  components  and  with  relatively  little  effort 
in  optimizing  and  fine  tuning  the  system  for  improved 
accuracy,  thereby  confirming  the  fact  that  accurate 
global  computation  can  be  performed  with  relatively 
inaccurate  individual  components. 

V.  Discussion 

The  number  M  of  state  vectors  of  length  N  that  can 
be  stored  at  any  time  in  the  interconnection  matrix  T y 
is  limited  to  a  fraction  of  N.  An  estimate  of  M  sse  0.1N 
is  indicated  in  simulations  involving  a  hundred  neurons 
or  less1  and  a  theoretical  estimate  of  M  =:  N/4  InN  has 
recently  been  obtained.2  It  is  worthwhile  to  consider 
the  number  of  bits  that  can  be  stored  per  interconnec¬ 
tion  or  per  neuron.  The  number  of  pixels  required  to 
form  the  interconnection  matrix  is  N2.  Since  such  a  T,; 
memory  matrix  can  store  up  to  M  »  N/4  InJV  (N- 
tuples),  the  number  of  bits  stored  is  MN  =  N2/ 4  In N. 
The  number  of  bits  stored  per  memory  matrix  element 
or  interconnection  is  M N/N2  =  (4  ln/V)-1,  while  the 
number  of  bits  stored  per  neuron  is  M N/N  =  M. 

The  number  of  stored  memories  that  can  be  searched 
for  a  given  initializing  input  can  be  increased  by  using 
a  dynamic  memory  mask  that  is  rapidly  addressed  with 
different  T y  matrices  each  corresponding  to  different 
sets  of  M  vectors.  The  advantage  of  programmable 
SLMs  for  realizing  this  goal  are  evident.  For  example, 
the  Litton  Lightmod  (magnetooptic  light  modulator), 
which  has  nonvolatile  storage  capability  and  can  provide 
high  frame  rates,  could  be  used.  A  frame  rate  of  60  Hz 
is  presently  specified  for  commercially  available  units 
of  128  X  128  pixels  which  are  serially  addressed.7  Units 
with  256  X  256  pixels  are  also  likely  to  be  available  in  the 
near  future  with  the  same  frame  rate  capability.  As¬ 
suming  a  memory  mask  is  realized  with  a  Litton 
Lightmod  of  256  X  256  pixels  we  have  N  =  256,  M  =* 
0. 1 N  se  26  and  a  total  of  26  X  60  =  1560  vectors  can  be 
searched  or  compared  per  second  against  an  initializing 
input  vector.  Speeding  up  the  frame  rate  of  the  Litton 
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Lightmod  to  increase  memory  throughput  beyond  the 
above  value  by  implementing  parallel  addressing 
schemes  is  also  possible.  Calculations  show  that  the 
maximum  frame  rate  possible  for  the  device  operating 
in  reflection  mode  with  its  drive  lines  heat  sunk  is  10 
kHz.7  This  means  the  memory  throughput  estimated 
above  can  be  increased  to  search  2.6  X  105  vectors/sec, 
each  being  256  bits  long,  or  a  total  of  6.7  X  108  bits/sec. 
This  is  certainly  a  respectable  figure,  specially  when  we 
consider  the  error  correcting  capability  and  the  asso¬ 
ciative  addressing  mode  of  the  Hopfield  model;  i.e., 
useful  computation  is  performed  in  addition  to  memory 
addressing. 

The  findings  presented  here  show  that  the  Hopfield 
model  for  neural  networks  and  other  similar  models  for 
content  addressable  and  associative  memory  fit  well  the 
attributes  of  optics,  namely,  parallel  processing  and 
massive  interconnection  capabilities.  These  capabili¬ 
ties  allow  optical  implementation  of  large  neural  net¬ 
works  based  on  the  model.  The  availability  of  nonlin¬ 
ear  or  bistable  optical  light  amplifiers  with  internal 
feedback,  optical  bistability  devices,  and  nonvolatile 
high  speed  spatial  light  modulators  could  greatly  sim¬ 
plify  the  construction  of  optical  CAMs  and  result  in 
compact  modules  that  can  be  readily  interconnected  to 
perform  more  general  computation  than  nearest- 
neighbor  search.  Such  systems  can  find  use  in  future 
generation  computers,  artificial  intelligence,  and  ma¬ 
chine  vision. 
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